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Abstract. Wc consider non-interacting multi-qubit systems as controllable 
probes of an environment of defects/impurities modelled as a composite spin- 
boson environment. The spin-boson environment consists of a small number of 
quantum-coherent two-level fluctuators (TLFs) damped by independent bosonic 
baths. A master equation of the Lindblad form is derived for the probe- 
plus-TLF system. Wc discuss how correlation measurements in the probe 
system encode information about the environment structure and could be 
exploited to efficiently discriminate between different experimental preparation 
techniques, with particular focus on the quantum correlations (entanglement) 
that build up in the probe as a result of the TLF-mediated interaction. We 
also investigate the harmful effects of the composite spin-boson environment on 
initially prepared entangled bipartite qubit states of the probe and on entangling 
gate operations. Our results offer insights in the area of quantum computation 
using superconducting devices, where defects/impurities are believed to be a 
major source of decoherence. 
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1. Introduction 

Superconducting qubits [1] consist of electronic nanocircuits embedding Josephson 
junctions whose dynamics can, in certain parameter regimes, be restricted to a two- 
dimensional manifold. These qubits can be used as test-beds for studying quantum 
mechanics at its most fundamental level [5], and are also potential candidates for the 
practical implementation of quantum information processors [3] . 

An interesting and challenging aspect of these endeavours is the process of 
decoherence, whereby qubits lose their quantum-mechanical nature and are rendered 
dynamically equivalent to their classical two-level counterparts. In the case of a 
superconducting charge qubit, such as the Cooper-pair box (CPB), dephasing (phase 
decoherence) is dominated by low-frequency noise thought to be caused by interactions 
with two-level fluctuators [H |5j [6l |7] (TLFs) in the local environment. These TLFs 
may be charge traps caused by defects/impurities in the Josephson junction or in the 
substrate. There has been a substantial amount of research on TLFs causing single 
qubit decoherence in Josephson junction systems. Theoretical works have concerned 
larger ensembles of TLFs, both incoherent [SlHIinillllinilTnilTl] and coherent [5], that 
randomly switch between two configurations, producing low-frequency fluctuations in 
the relevant qubit parameters. As discussed in detail in all these works, whenever 
a large number of TLFs are very weakly coupled to the qubit, their effect can be 
described by a conventional boson bath with a suitable chosen spectral density. This 
is however not always the case — there may be situations when only one or a few 
impurities are important. Indeed, Neeley, et al [H] have demonstrated the existence of 
coherent TLFs. Lupa§cu, et al ^ provide evidence that these TLFs are in fact genuine 
two-level systems. In Neeley's experiment [B] a single TLF that was coupled to the 
qubit led to an avoided crossing in the qubit energy spectrum. The TLF was used as 
a proof-of-principle memory qubit, but such TLFs will in general be detrimental to 
the operation of superconducting qubits. It is therefore desirable to understand the 
behaviour of TLFs in the vicinity of superconducting qubits in order to better-equip 
quantum information scientists to manage the challenge of decoherence. 

In view of its importance for understanding decoherence in superconducting 
nanocircuits we study in this paper the effect that a few underdamped, coherent 
TLFs have on the quantum dynamics of a charge qubit. A single qubit coupled 
to a single coherent TLF, that was in turn (under-) damped by a bosonic bath was 
already considered in [15i . We find this type of non-equilibrium environment [16j 
to be an appealing model of the environment of a Josephson junction qubit because 
TLFs may be damped by phonons in the substrate, for example. We refer to this 
as a composite spin-boson environment. The case of a few TLFs, which will be 
considered here, is rather complex and therefore we resort to extensive numerical 
calculations. Our approach represents quite a general strategy for probing a composite 
spin-boson environment, whereby we derive a Markovian master equation for the 
average dynamics of the probe-|-TLFs after tracing out the external baths. To allow 
for the random nature of TLF formation, we select the TLF Hamiltonian parameters 
according to probability distributions designated in [5J. Of the many environmental 
properties that could be studied in this scenario, we focus on inferring the presence 
or absence of coherent coupling between the TLFs. This connectivity of the TLFs has 
been identified as important in the decoherence of a bipartite qubit system [17 , with 
further-reaching importance for quantum computing. 

Performing a more thorough treatment of this situation is increasingly important 
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in light of recent experiments [6j |7j involving superconducting qubits coupling 
predominantly to only a small number of TLFs. In these experiments, measurements 
on the qubit were used to infer properties of the TLF (although that wasn't the focus 
in That is, the qubit was used as a probe of the environment |THl|Tn]. Probing 
properties of a small collection of TLFs, such as whether or not they interact with 
each other, necessitates consideration of the full dynamics of the probe plus TLFs 
under the influence of external baths. 

From the theory point of view, identifying structural properties in the 
environment could be done using some form of noise correlation measurements. Here 
we focus on inferring environmental features via the analysis of the entanglement 
that will build up in a probe consisting of two non- (directly) interacting qubits whose 
remote coupling is mediated by the TLFs present in the surroundings. When the 
probing is "local", so that each probe qubit couples to just one single TLF, the 
absence of quantum correlation generation would immediately signal a non-connected 
environment, given that the probe qubits can only become entangled if the fluctators 
would couple to each other. Entanglement swapping in those circumstances has 
been discussed in the literature [201 EH [22] • When probe qubits are subject to 
the action of a few TLFs, as it happens in qubit realizations in the solid state, 
we will show that the remote entanglement in the probe bears signatures that 
can be linked to the connectivity in the environment and can in some cases be 
related to monogamy constraints Given that bipartite entanglement has been 

shown [21] to be lower-bounded by combinations of pseudo-spin observables, we also 
analyze what information can be extracted from magnetization measurements along 
a given direction (in the case considered here the magnetization along the z-direction 
corresponds to the average charge) and study the power spectra of magnetization 
observables using both single and bipartite probes. We find that a double-qubit probe 
generally outperforms a single-qubit probe, a result that could perhaps be expected 
given the extra degrees of freedom available in the composite system. We supplement 
our analysis of correlation measurements by investigating the decoherence of composite 
probes initially prepared in a certain maximally entangled state when subject to a 
composite spin-boson environment, as well as the performance of entangling gate 
operations when performed in the presence of this type of noise. 

The paper is organized as follows. Section [2] sets the scene by describing our model 
for the double-qubit probe and spin-boson environment. A detailed derivation of the 
proposed master equation as well as a discussion of its validity domain are presented 
in the appendix. Numerical results for probing the connectivity of the spin-boson 
environment are presented in section [3] using both probe entanglement and estimated 
power-spectrum analyses. We summarize and discuss these results in sectio n [3^ 
as well as compare the double-qubit probe to a single-qubit probe. In section |4] we 
investigate the decoherence of maximally-entangled Bell states induced by a composite 
spin-boson environment. The performance of bipartite entangling gates in the presence 
of this form of noise is analyzed and discussed in section [s] Section [6] concludes. 

2. System 

The system we consider is illustrated in figure [T] It consists of two charge qubits 
(blue spheres) acting as probes of an environment containing TLFs (grey spheres). 
Each qubit is coupled to a few TLFs (black lines in the figure) but probe qubits are 
assumed to not directly couple to each other. In the numerical calculations we will 
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consider the case in which there are four TLFs. Four TLFs is a balance between 
generating the desired spectral features (requiring an ensemble of TLFs [31 [H] ) and 
maintaining reasonable computation time (smaller Hilbert space). Also, it may be the 
case that only a few TLFs will couple strongly to a Josephson junction qubit, as in 
recent experiments |BJ [7] . We should stress that the conclusions of our work do not 
depend on this choice. 

In the charge basis, each qubit/TLF has a local free Hamiltonian consisting of 
both longitudinal (a^) and transverse (tTx) components: 2H„ — ea^ + Actx. In the 
eigenbasis the corresponding pseudo-spin Hamiltonians are 2Hs — ^aSz, where the 
spin frequency is ^1 — + A^. (Throughout this article we denote Pauli operators in 
the charge basis by tr, and in the pseudo-spin basis by s.) For simplicity, we "engineer" 
the probe qubit Hamiltonians to have only longitudinal components (Ap — 0). 

Probe: We label the identical probe qubits A and B. Choosing uncoupled probe 
qubits for reasons that will become clear later, the total Hamiltonian for the probe is 
the sum 2Hp = rip(s^ + s^)- 

Impurities: We label the four TLF impurities with j = 1,2,3,4. The total 
Hamiltonian for the TLFs is then 2ffTLF = J2'j=i f^^Sz"''' + Vtlf, where Vtlf describes 
coherent couplings between the TLFs, if they exist (defined below). Note that the 
pseudo-spin basis for the impurities is different to the probe (the Sz'^ axis is rotated 
relative to Sz) because the probe and TLF energies will differ, in general. Recent 
theoretical work [5] suggested specific distributions of these TLF energies in order to 
account for both low- and high-frequency noise observed in superconducting quantum 
systems. In our numerical study, we have adopted these distributions to determine 
the TLF bias energies ej (linear distribution) and tunnel amplitudes Aj (log-uniform 
distribution) . Throughout the paper we will refer to Aj as the local field. Again, note 
that our results are independent of the specific choice of frequency distribution and 
the same qualitative results can be derived when using a different functional form, e.g., 
a linear or a uniform distribution in a selected interval around the qubit frequency. 



Figure 1. Double-qubit probe schematic. Blue splieres are tiie probe qubits, 
Alice and Bob. Grey spheres are the TLFs. Probe— TLF coupUngs (uj) are 
depicted as black lines. TLF-TLF couplings {fJ.j^k) a-re specified by red lines. 
Each fiuctuator is also subject to the action of a bosonic bath at a temperature 
T. Interactions between A and B are mediated by the TLFs, which lead to remote 
(TLF mediated) entanglement generation in the probe. 
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Interactions: It is sensible to expect that the dominant interaction in a system 
of coherent two-level charges is an electrostatic one [101 US]- That is, charge- 
charge interactions. We therefore assume bipartite ZZ interactions (ctz (E) CTz) between 
subsystems. Within the TLFs we assume nearest-neighbour ZZ interactions of 
strength /ij.fe, where k — {j mod 4) + 1. We assume no direct interaction between 

the probe qubits A and B. So, Vtlf — Mj,fc'5'z''^'^z'^^- For coupling strength 

between the jth impurity and the probe qubits, we have Vp = Sj=i ^A^^ + Sr^)ai^\ 
We define V = Vp-I-Vtlf- In the numerical simulations, noting that we expect distant 
TLFs to have very little impact, we have assumed all couplings ^.j^k = fi and all Vj = ly 
to avoid unnecessarily cumbersome results. Our conclusions are valid even when small 
variations in the parameters, of the order of 5 — 10%, are considered. 



2.1. Master equation 

The impurities are coupled to independent reservoirs of bosons (e.g., phonons in the 
substrate), leading to dissipation (damping) as in the spin-boson model [26 . Under 



appropriate weak-coupling assumptions (see Appendix A for our derivation), the 
dissipative dynamics of the composite system (qubits plus damped TLFs) can be 
expressed in the following Born-Markov master equation for the joint state p{t) of the 
probe plus impurities: 

m = -i[H, Pit)] + ^(Pp) + + I?i^'V(i), (1) 

j 

where the total Hamiltonian is 

H^Hp + Htlf + V. (2) 

The V superoperators |27j represent decoherence in the TLFs due to coupling with 
the bosonic baths, which are at temperature T. The decoherence consists of dephasing 
Vi'^p = Ti'\si^'^ psi^^ - p], emission into the baths Pi^'V = T^-^s'-l^ ps^+ - (s^'^s^ V + 
/3S^^s^'V2], and absorption from the baths v'^^ p = F^j!'' [s^'^s^' - (s^^s+V + 
ps^^s^'')/2]. The TLF decoherence rates r^"'^- are proportional to the respective 



dephasing, emission and absorption rates 7z.± (see Appendix A), which are functions 
of the temperature and spectral properties of the jth bosonic bath (with absorption 
dramatically reduced at low temperatures). Further, the TLF decoherence rates are 
also functions of the ratio of local field to bias, tanOj = ^j/sj [28 . To get a feel 
for the influence of this ratio, if we assume dissipation- limited dephasing (72 — 7-/2) 
and sufficiently low temperature (7-1- /7_ 0), then tan^j dictates the dominance 
of pure dephasing or relaxation in each TLF. Specifically, ri"'''/r*f '' = l/tan^6'j 
so that pure dephasing dominates the TLF decoherence for weak local fields, and 
relaxation dominates for strong local fields. Following |5;, we distribute the random 
TLF parameters Sj, Aj, and 72'^j. as per the distributions P{ej) oc Sj, -P(Aj) cx 1/Aj, 
P(72"'j.) c>c I/7P].. We choose these parameters to take values within the following 

moderate ranges: Ej e (l±0.5)ej; Aj e Aj±0.5 min(rip, Aj); 7^''^. G [rimin/S, rimin/2], 
where f^min is the minimum spin frequency amongst the TLFs. We are free to select 
sensible values for the overbar quantities Ej and Aj, which we will reference to the 
tunable probe frequency Hp. Importantly, the TLFs are underdamped (therefore 
requiring a quantum- mechanical description), so 7'-'-* < flj. We take the probe- 
TLF coupling to be uniform (vj = v) and weak compared with all of the TLF 
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frequencies: v — rimin/S. Our assumption of weak probe-TLF coupling simplifies 
the master equation derivation significantly (see equation A.l a detailed discussion 
of an analogous situation can be found in [H]), and is in addition to the standard 
Born-Markov approximation of weak TLF-bath coupling. Weak probe-TLF coupling 
is in accord with the recent experiments of [HI |7] , as well as the experiment of |1] (see 
|9]) where I'j ^ MHz and fij ^ fip ~ GHz. The TLF-TLF coupling is also assumed 
to be uniform ^j^k — Note that the authors of [5] point out that v must also be 
randomly distributed in order to realize 1// noise in the probe. Generating the correct 
statistics is not within the scope of this paper as we are interested in describing effects 
when the environment is dominated by only a few fluctuators. 



2.2. Observable quantities 

We restrict our knowledge to the probe subsystem (as would be the case in an 
experiment). A notable observable quantity on the probe is the magnetization, 
which is related to a simple sum of Pauli operators: M^{t) = + s^. The 
appeal of considering the probe magnetization is that it requires only tractable, 
local measurements on each probe qubit. That is, our results may be easily tested 
in an experiment. It is worth noting that a result of Audenaert and Plenio |24j 
shows that measuring correlations C^^j^^ = ((t^/z ® ^xjz) along the XX and TIL 
'directions' suffice to give a lower bound on the probe entanglement. This can remove 
the requirement for full tomographic (entanglement) measurements when verifying or 
quantifying entanglement in the probe. 

The time series resulting from measuring the probe's X-magnetization is M(t) = 
{M-^(t)). The mean-square power spectrum of M{t) is given by 

/oo 
i?(r)e--Mr, (3) 
-oo 

where the reduced auto-correlation function is R{t) = {M{t + T)M{t)) — {M{t H- 
T)){M{t)). The angle brackets here denote the expectation value of an operator 
{x) = Tr[xp(t)]. In reality, M{t) is a discrete quantity (the data is a time series), 
and so the power spectrum obtained is an estimate (the Fourier transform of the 
reduced autocorrelation of the time series M{t)). This estimate of average power 
as a function of frequency can theoretically be improved (by increasing the duration 
of the experiment, for example), but this may not be practical in reality. Strictly 
speaking, M (t) must be a wide-sense stationary process (time-independent first and 
second moments) for the power spectrum to exist. 



3. Results: Detecting the presence of coupling between the TLFs 

Can probe observables reveal the degree of connectivity of a composite spin-boson 
environment? In this section we present numerical results showing that this is 
indeed the case. Observable quantities we consider are the probe magnetization, its 
estimated power spectrum, and the remote entanglement between the probe qubits. 
Entanglement generated between probe qubits that are initially in a separable state 
is primarily due to the structure of the spin-boson environment (e.g., the presence or 
absence of TLF-TLF interactions in the surroundings of the probe). 

Initially we set the probe to be in a state orthogonal to the eigen-axis: \ipp{0)) — 
|+)a|+)b- The TLFs are assumed to be initially in a zero-temperature thermal state. 
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i.e., the ground state of Htlf, which we denote as \g). So, |0tlf(O)) — \g)- We 
plot observable quantities as a function of the ratio of TLF-TLF coupling strength to 
probe-TLF coupling strength, fi/i/, which is often believed to be small (see |TH], for 
example). 



Section 



3.1 



considers the estimated power spectrum of (M^it)). Section 



3.2 



considers the build up of entanglement in the probe. 



3.1. Power spectrum of (A/x) 

3.1.1. TLFs with weak local fields. Consider the case of "weak" local fields in the 
spin-boson environment where tanOj ^1/3 (specifically Aj/sj = 1/3, where TLF 
dephasing dominates relaxation: Fz/F- ~ 10). In our probing of TLF connectivity, 
we tune the ratio of probe-TLF splitting ep/ej to 10, 3 and 1, corresponding to 
figures p(a)| 2(b)| and 2(c)[ respectively. We observe two effects: 1. For weaker TLF 



interconnectivity ii/v < 0.6, a decrease is observed in the height of the single dominant 
peak in the spectrum; 2. For stronger TLF interconnectivity ii/v ~ 1, the spectrum 
splits into multiple peaks, the most dominant of which is shifted in frequency relative 
to ii/v = 0. This is visible in figure [2] where the gold traces (ii/i' — 1) are the most 
qualitatively different from the blue traces {fi — 0). The power in the signal {M^{t)) 
redistributes from one dominant frequency for unconnected TLFs (fi — 0), to multiple 
frequencies as ^/v approaches unity (highly connected TLFs). Remarkably, the most 
dominant peak for highly connected TLFs ji/v — \ is qualitatively similar to the case 
of an isolated probe (figure |2(d)| — a single peak a,t uj — ilp — although the peak 
visibility (height) is noticeably less than the isolated probe. 



3.1.2. Effect of TLF local field. Above we have established that we can distinguish 
between highly connected TLFs (/i = v) and unconnected TLFs (/i — 0) for weak 
local field strength Aj < Sj. It is interesting to ask how stronger local fields affect 
our ability to distinguish these two values of /i. To explore this, we increase the ratio 
of TLF field strength to splitting: tan^j — Aj/ej, which was less than 1 in figure [2] 
This varies the range for Aj, which (we remind the reader) we have taken to be 
Aj ± 0.5 min(Aj, rip). This range ensures a sensible variation in the TLF local field 
strengths of no greater than one-half of the probe frequency. (We assume that the 
TLF energies are distributed over a relatively small range as might be expected for 
systematically formed impurities/defects.) 

As the local field increases tan^j > 1, two different effects occur: 1. Larger local 
fields cause relaxation to dominate over dephasing in the TLF decoherence; 2. The 
TLF eigenstates increasingly align towards the axis, and seem to have a decreasing 
effect on the probe, perhaps because the interaction is of the CTz (E" CTz type. Evidence 
to support this is shown in figure [s] which shows a peak visibility (height) reduction 
of about one order of magnitude as tan 9j is increased by one order of magnitude from 
1/3 to 3 |3(a)| to [3(b)| . So, although extra features appear in the power spectrum, they 
become increasingly difficult to observe. Despite this, there is at least one plot in each 
column for which /i = and fj, = v are distinguishable. Thus it is apparent that tuning 
the probe frequency (selecting a row in figure [3]) allows these TLF connectivities to be 
distinguished for a wide range of values of the TLF local field strength (we obtained 
similar results for tan^j = 1). 
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(d) Control spectrum: u = 0. 



Figure 2. Estimated power spectrum of {M^it)) for weak local fields tanSj = 
1/3 with the initial probe state |'i/'(0)) = |H — h). Sa mplin g parameters (units of 
rip = ep): (5/ = 1/200, is = 0.05, /jv = 10. Figure |2(d)] is the power spectrum 
for an isolated probe. The insets show the smaller peaks magnified. The relative 
heights and positions of the peaks can be used to distinguish unconnected TLFs 
(/.t = 0, thick solid blue line) from highly connected TLFs (/i — > v, thick dashed 
gold line). 



3.2. Probe entanglement 

Interactions between probe qubits are mediated by the TLFs and entanglement 
between the probe qubits can be generated in this indirect way. We now consider using 
the probe entanglement to distinguish between connected and unconnected TLFs. We 
use the logarithmic negativity as a measure of bipartite entanglement between the 
probe qubits, defined as [3T] 

£;p=log2||p^-||i, (4) 

where || • ||i denotes the trace norm, and p^* is the partial transpose of p. 

3.2.1. TLFs with weak local fields. Figure |4] shows the logarithmic negativity of 
/9p(i) for the same data sets as in figure |2] (weak local fields tan^j = 1/3). We make 
two observations. Firstly, tuning the probe frequency provides only one benefit for 
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Figure 3. Estimated power spectrum of {Mx(t)) for ep/ej = 10 (top row), 
3 (second row ) and 1 (t hird row). As the local field strength tandj = ^j/sj 



increases from 
discussion 



3(a)(]3(b)| TLF relaxation dominates TLF dephasing. See text for 



distinguishing /i = from /i 7^ 0. This is evidenced by the remarkable quahtative 

4(c)[ The benefit is the number of probe 
qubit cycles required to distinguish the two cases (time is shown in units of probe 
qubit cycles). Secondly, it is clear that strongly coupled TLFs (the gold line) cause 
less entanglement to generate within the probe. This can be explained by invoking 
the concept of entanglement monogamy [531 ISS]- As the TLF-TLF connectivity 
increases, the indirect link between the two probe qubits is weakened, and so remote 
entanglement generation slows. It is important to emphasize which partitions one 
should consider when invoking monogamy arguments. The relevant quantity is 
the entanglement shared between each probe qubit and a given TLF; this is the 



similarity between figures |4(^ gbM and 
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quantity that is sensitive to entanglement sharing in two ways: (1) It decreases as 
the entanglement in the probe builds up, independently of the connectivity in the 
TLF environment; (2) It "feels" the TLF-TLF coupling in the sense that, within a 
selected time interval, the stronger the fluctuators couple, the smaller the entanglement 
between probe qubit and TLF becomes. As a result, remote entanglement builds 
up more slowly when fluctuators couple so that it quickly degrades in a decohering 
environment, as illustrated by results in figure [5] 
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(a) £p = lOsj. 



(b) ep = 3ej 



(c) ep 



Figure 4. Entanglement (logarithmic negativity) in the probe as a function of 
probe cycles, for small TLF local field strength tanSj = 1/3. Note the distinct 
diff'erence between strongly interacting TLFs (fi u) and non-interacting TLFs 

(m = o). 



3.2.2. Effect of the TLF local field. Unlike the power spectrum (figure [s]), the probe 
entanglement Ep{t) remains useful for distinguishing /i = from fJ. ^ for strong 
local fields in the TLFs where tan^_,- > 1. Figure [s] shows the probe entanglement as 
a function of time for the same data sets as in figure [3] It is clear that entanglement 
is generated within the probe for some time, before the TLF decoherence causes it to 
dissipate. As argued before, this loss of generated probe entanglement occurs faster 
for highly connected TLFs with fj./v = 1, as one might expect, even without no 
explicit mention of monogamy constraints, since these TLF-TLF connections provide 
more links between the probe qubits and the TLF decoherence channels. This 
faster dissipation of generated entanglement for highly-connected TLFs allows us to 
distinguish between fi/v = Q and ji/v ^ \ after 10 to 50 probe qubit cycles, depending 
on the TLF parameters. Similar conclusions can be drawn for tan^j = 1. 

3.3. Discussion of the results 

Since we are performing local measurements on each probe qubit, one might ask if 
there are any advantages of using a double-qubit probe. An obvious advantage is that 
entanglement within the probe becomes an accessible quantity that is not possible in 
a single-qubit probe. This is important for detecting the TLF-TLF connectivity, as we 
have seen that this task was achievable over a wider range of TLF parameters using 
the probe entanglement than the power spectra of the probe magnetization (which 
yielded qualitatively similar results for both types of probes). 

In section [3] we were able to distinguish between composite spin-boson 
environments with high connectivity from those with low connectivity. The parameter 
ranges for which our findings were robust are: ep/ffj = 1,3,10 (tunable), 1/3 < 
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Figure 5. Entanglement (logarithmic negativity) between the probe qubits for 
Sp/sj = 10 (top row), 3 (second row) and 1 (third row). The cases fi = (solid 
blue line) and fi = u (dashed gold line) are qualitatively distinguishable unless 
both Ej/ep = 1 and tanSj > 1. Note that the upper limits of the axes change 
between plots. Same data sets as in figure Is] 



t'dnOj < 3. Other parameters (F, etc) are restricted to lie within the ranges that 
ensure validity of the master equation — see [Appendix A| 

For all numerical calculations in this paper we have assumed effectively zero- 
temperature bosonic baths where fij — [exp(hflj /kBT) — 1. The relevant 
frequencies for experiments with Josephson qubits are in the vicinity of 10 GHz 
[H |32l ini |7], with cryostat temperatures of the order of 30 mK 2]. These values 
give n 0.1, so the low-temperature approximation is good. 

In general, determining the probe entanglement would require full quantum-state 
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50 100 150 200 50 100 150 200 

(a) 11 = 0. (b) n = V. 



Figure 6. Lower bounds on the probe entanglement as described in 24 . Only 
the optimal bound (dotted red line) provides a good approximation to the 
probe entanglement. Parameters are tanffj = 1/3, ep/sj = 1, and |V'(0)) = |H — h)- 



tomography. Here we consider estimating the probe entanglement from measurements 
less costly than full quantum state tomography, as described in (and references 
within) . The result is lower bounds on the entanglement — we refer the reader to |24j 
for details. 

For all parameter regimes considered in this paper we found that one of the lower 
bounds given in [24 — the optimal one given in (|5| — provided a remarkably good 
approximation to the probe entanglement for all times. The other lower bounds Ci^2 
in 24J didn't approximate the probe entanglement well for any time. An example is 
shown in figure [6j The solid line is the entanglement within the probe (logarithmic 
negativity), and the dotted line shows the lower bound given by 

C^(pp) = max[0,log2(l + |Ai| + IA2I + IA3I) - 1], (5) 

where Ai^2.3 are the eigenvalues of the matrix 

/ c"="= c^y c^^ \ 
A = 0^" . (6) 

\ c^^ c^f c^^ / 

The matrix A is formed from probe observables c*^ = Tr[sf^ ® ^^Pp] ~ 3;,y,z). 
Note that c*^ — c'* due to the symmetry of the problem. 

4. Decoherence of entangled states 

Entanglement has been identified as a key resource for quantum information 
processing. It is therefore important to study the loss of entanglement induced by 
coupling with an environment, as this coupling is generally unavoidable. In this section 
we reconsider our double-qubit probe as a double-qubit register (DQR) interacting 
with the same spin-boson environment (four damped TLFs) as above. Starting the 
DQR in an entangled pure-state, and the TLFs in their ground state (as before) 
with weak local fields tan^j < 1, we numerically investigate the behaviour of the 
logarithmic negativity as a function of time. Specifically, we consider the lifetime of 
distillable entanglement (for which the logarithmic negativity is an upper bound) . We 
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compare our results to previous studies of entanglement decay [Ml IM], all of which 
used rather less sophisticated models for the environment. Nevertheless, we find some 
qualitative similarities between our results and previous work. 

Reference |33j consider multiqubit states whereby each qubit is damped by 
independent baths. We refer to this as 'direct' damping, by an equilibrium 
environment (a reservoir). In our non-equilibrium spin-boson environment, the 
damping is mediated by the TLFs and we refer to this as 'indirect' damping of 
the DQR. Reference [33] parameterizes time via the probability for a qubit to 
exchange a quantum of energy with its bath (in the absence of pure dephasing), 
p(t) = 1 — exp[— 7(2n + l)t/2\. Here 7 is the zero-temperature damping rate, and n 
is the mean number of excitations in the bath (n is zero temperature). When 
the DQR logarithmic negativity falls below an arbitrarily small fraction of its initial 
value, e <C 1, the distillable entanglement can be considered zero. The time at which 
this occurs is t^. For generalized GHZ states (requiring > 3 qubits), and for three 
different types of direct damping, [57 found that p{te) oc — log e (although they were 
looking at it from a slightly different perspective, as we comment later). Remarkably, 
we find numerical evidence of the same qualitative behaviour for the decay of the 
Bell states in a DQR (see figure [7|. 

For low temperature T ~mK (as we have considered throughout this paper unless 
otherwise noted), weak TLF local field tandj = 1/3, and TLFs with relatively small 
charge: sp = 3ej, we consider the DQR to be initially in each of the four Bell states 
in turn: |</.±) = (|00) ± \n))/V2, \^^) = (|01) ± |10))/V2. Figure [t] shows the DQR 
logarithmic negativity as a function of time, as well as p{te) = 1 — exp(— te/2) for a 
DQR initially in the states It is evident that p{t^) oc — log e for \<f>^). Further, we 

can see that interacting TLFs (the gold traces) tend to reduce the DQR entanglement 
faster. This is expected since interacting TLFs provide more connections between the 
DQR and the baths. For the chosen set of parameters (particularly A = for the 
DQR), the states 1-0*) commute with Hp + Vp and so don't evolve, nor couple to 
the TLFs. We found the same qualitative behaviour for ep/ej = 1 and 10 when the 
TLF local fields were not strong: tandj < 1. For strong local fields tan^j > 3, the 
linear relationship p(<e) (x — log e did not hold in general. This was because the probe 
entanglement tended to exhibit quite erratic behaviour, such as multiple collapses and 
revivals. 

A comment on the previous work in |33l 134] is appropriate here. In those works, 
the decay of TV-particle entanglement was considered as a function of TV. In [33] 
it was found that p{t) cx — (l/A^)loge. Here we have fixed N = 2 and found that 
p{t) oc — loge. Our focus is slightly different, but it is interesting that the loss of 
distillable entanglement (logarithmic negativity) is qualitatively the same for direct 
and indirect damping of bipartite qubit states (within the parameter regimes discussed 
in the previous paragraph). It is important to remark that the coincidence with 
the predictions for the decoherence of multipartite states subject to independent 
reservoirs should not be considered as a general result given that we analyzed a 
very special case, which is the one of two entangled qubits in selected parameter 
regimes. What is relevant for our purposes is the fact that the agreement with the 
analytical prediction in |33j for direct decoherence points out a sharp asymmetry in 
the processes of entanglement "destruction" and (remote) entanglement generation 
in a composite environment, so that there are circumstances where the TLF systems 
may be essentially invisible when analyzing the decoherence of initially entangled 
probe states, while the presence of the TLFs would be revealed when monitoring 
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Figure 7. Decay of entanglement (logarithmic negativity) in a double-qubit 
register initially in the Bell states \<j)'^) (top row) and \4>~) (second row) for 
uncoupled TLFs {fi = 0, solid blue line) and coupled TLFs {fi = dashed gold 
line). Other parameters are tanOj = 1/3 and e-p/sj = 3. See text for discussion. 

entanglement creation in the probe. 

5. Effect of spin-boson environment on entangling gate operations 

In this section we investigate the effects of the composite spin-boson environment 
on the performance of entanghng gates. Starting the DQR in the separable state 
|'!/'p(0)) = we consider two entangling gates: a ZZ gate (s^ ® sf); and an 

XX+YY gate [s^ <S) + Sy Sy]. In the ideal case there are no TLFs and bipartite 
entanglement (quantified again by the logarithmic negativity) is generated between 
the isolated register qubits in an oscillatory fashion as shown by the solid red curves 
in figure |8] In the presence of four TLFs, the entangling gate performance is clearly 
reduced, and further modified depending on the strength of the local fields in the TLFs. 
This may be understood as follows. We have argued that the presence of a coherent 
coupling between the TLFs leads to a decrease in the effective interaction strength 
between the qubits in the probe, i.e., the "effectiveness" of the indirect link between 
probe qubits is diminished, a result that can be interpreted in terms of monogamy 
constraints leading to a slow down in the process of remote entanglement creation. 
In the case of the entangling gates, one could perhaps argue similarly, considering 
now a bipartition separating the probe qubits. The higher the connectivity in the 
environment, the slower the entangling gate can operate, as illustrated in figure [9] 
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(a) tanej=l/3. (b) tanSj = 3. 



Figure 8. Entanglement (logarithmic negativity) between the register qubits for 
\ip{0)) = |++>, ep/gj = 1 generated by XX+YY (top row), and ZZ (bottom row) 
gates. The presence of TLFs (blue dashed and green dotted lines) diminishes the 
performance of the entangling gates, as would be expected. 

At longer times (the order of 100 register-qubit cycles and greater), the DQR 
appears to approach an entangled steady state for weak local fields tandj = 1/3. This 
is shown in figure |9]for the ZZ gate (the same qualitative behaviour occurred for the 
XX+YY gate). We were unable to obtain analytical results to verify the presence 
of steady-state entanglement in the DQR generated by either gate (ZZ or XX+YY), 
when coupled to four TLFs. Our numerical study found that entanglement generated 
in the DQR was dissipated more rapidly (in less probe qubit cycles) for strong local 
fields in the TLFs. 

6. Conclusion 

We have considered superconducting qubits which are subject to decoherence 
dominated by low-frequency noise thought to be produced by interactions with a 
small number of defects/impurities. We model these impurities as coherent two-level 
fiuctuators (TLFs) that are under-damped by baths of bosonic modes (e.g. phonons). 
We probed such a composite spin-boson environment by making measurements on 
a pair of noninteracting qubits that each interact directly with the TLFs. Our 
extensive numerical study revealed that the presence or absence of coherent coupling 
(connectivity) between the TLFs can be discriminated in two ways: from the 
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Figure 9. Long-time entanglement (logarithmic negativity) between tlie register 
qubits generated by a ZZ gate, for |</'{0)) = |H — h), sp/Sj = 1 and tanfij = 1/3. 

estimated power spectrum of the probe magnetization (requiring relatively long-time 
measurements) and from entanglement generated within the probe, mediated by the 
spin-boson environment. We argue that entanglement monogamy considerations |23j 
allow interpretation of our results in terms of an effective decrease in the remote 
interaction strength when the environment is connected, which yields to a creation of 
quantum correlations on a much larger time scale as compared with the uncoupled 
fluctuator case. 

We also showed that this remotely generated entanglement can be well-estimated 
by a lower bound [24] that requires less experimental effort than the full quantum 
state tomography required to evaluate the entanglement. The upshot is that this 
connectivity of the TLFs should be discernible using tractable measurements in a real 
experiment. This result is important for studies of quantum- mechanical phenomena 
in Josephson devices (including quantum computing) where it is desirable to minimize 
the effects of decoherence, for which the TLF-TLF connectivity may play a significant 
role [T7]. 

When considering the effects of the spin-boson environment on a double-qubit 
register initially prepared in a maximally-entangled (Bell) state, we showed that the 
presence of TLFs may be unnoticeable in certain parameter regimes, in the sense that 
entanglement degradation there is well approximated by the same decrease law as for 
direct decoherence. This fact emphasizes the possible usefulness of monitoring the 
reverse process of entanglement generation for environmental probing. 

The presence of TLF-TLF coupling also showed a reduction in the performance 
of entangling gate operations performed on the register. Our extensive numerical 
study can be supplemented by an analytical treatment of a simpler situation valid for 
short times, where TLF decoherence can be ignored. These results will be presented 
elsewhere [55] . 
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Appendix A. Derivation of the master equation 

Consider a single charge qubit "system" coupled to a finite number of independent 
charged impurities that fluctuate coherently between two configurations. These 
coherent two-level fluctuators (TLFs) are coupled to independent bosonic baths 
(phonons in the substrate, for example) that produce damping. The total Hamiltonian 
is the sum of free Hamiltonians for the single qubit, the TLFs, and the baths, and the 
interaction Hamiltonians: 

H = Hsq + HtlF + Hq + T4q-TLF + ^^TLF-B- 

The free Hamiltonians are H^q = (e^z + Aa^)/2, Htlf = Y.ji'^j^^^''' + 

Hb = i^i'^i jf^i j (the baths can be assumed to be identical so uji,j = uig). 

The coupling Hamiltonians are t^q-TLF = '^ji^jO'z'''o'z (Coulomb interactions), 

Vtlf-b = i ^iO'z'\o-g j + al j) (similarly to LOi, the couplings are assumed to be 
independent of the baths: Xe.j = Xe)- We denote the charge-basis Pauli operators by 
o'x.y.z- At this stage the TLFs are not interacting. 

The TLF-related energies Sj, Aj and Vj are randomly distributed following 
independent distributions discussed in |5] . The details of these distributions are critical 
for realizing the experimentally observed 1// noise spectrum of the qubit voltage/bias. 

We refer to the eigenbases of Hgq and -ffxLF as the pseudo-spin bases, and denote 
the corresponding Pauli operators as Sx,y,z- The Hamiltonians in the pseudo-spin basis 
are 



where tanSj = Aj/ej. In order to obtain the Master equation we move to an 
interaction picture with respect to the Hamiltonian 

Ho = + Htlf + Hb + Kq-TLF- 
In this picture the evolution equation for the total system is: 



Htlf 



H, 



sq 




i 




dptotjt) 
dt 



i[VTLF-B(i),Ptot(t)] 
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where 

Ptot(i) = exp(zJ?o<)ptotexp(-iJ?o^), 
VTLF-B(i) = exp(iJ?o<)VTLF-Bexp(-i7?o^)- 
By iterating once the above equation as usual [Ml [55] we obtain 

dm 



dt 



= -Tre„v / di'[V^TLF-B(i),[V^TLF-B(i'),Ptot(t')]] 

Jo 



where Tronv denote the trace over all of the baths. By assuming a factorized initial 
state of the form /9tot(0) = p(O)0j Pfc.ji with p(0) the initial state of the "qubit + 
TLFs" and pbj a thermal state of the jth bath, we may make a Born approximation 
in the coupling constants A^. The evolution equation becomes 

^ ~ -^^on.J^ di'[l>TLF-B(i),[VTLF-B(i'),P(i')(8)/^M]]- 

Since the baths are independent and they are all in a thermal state (diagonal in the 
number basis), by inserting the expression for VTLF-B(i) it is easy to check that 



Tr_V tdt' 



dt Jo 



\i,{cosej~si^\t) - sinefs^^\t)){hi,^j + ~a\^), 



3 

This is basically a sum of the expressions obtained in the standard derivation [3H1 [23] 
for the one qubit case. Now we make the next crucial assumption in the derivation. 
Assuming weak qubit-fluctuator coupling whereby Vj <C (as in the experiments of 
[H [SI [7] — see below) , we may approximate the Heisenberg operators of the TLFs in 
the interaction picture as: 

s{t) ~ exp(«/fooi)sexp(— i//oo^)i (A^-1) 

where Hqq = Hgq + H^lf is the sum of the free Hamiltonians. For example, we find 

that s^l\t) = s^£e^'^^K 

One may question the validity of our assumption of weak qubit-fluctuator 
coupling. For guidance we consider the experiments of |1] (charge qubit), [5] (phase 
qubit) and [7 (flux qubit), where i^j < O.lGHz and ^Ij ^ U ^ GHz, which falls within 
this weak-coupling regime (the qubit-fluctuator coupling strength in the experiment 
of ^ is estimated in (^). In this regime we can insert the result of the one qubit case 
for every fluctuator j: 

+ rL^) {s^l\t)~p{t)s^l\t) ~ [s^i\t)s^l\t)p{t) + ~pitr4\t)'^^-\t)]/2 

+ T^l^ {S^i\t)p{t)s'--l\t) ~ [s(f)(t)stHi)p(<) + P(i)s(f'^(i)sf Hi)]/2 

where the TLF ladder operators are s± = (sx ± iSy)/2, and the decoherence rates 
are V^^ = -i^^cos^ 9^ /2, r*f'^ = (7_ -I- 7+ ) sin^ 6*^ /4 and F^'^ = -i^sin^Oj/A. Here 
7z, 7_ and 7_|_ are the dephasing rate, the spontaneous emission rate and stimulated 
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emission rate respectively, which can be calculated by knowing the spectral properties 
and temperature of the bath. 

Now we return to the Schrodinger picture. For consistency the Heisenberg 
operators need to be sent back into the Schrodinger picture with the same -ffoo used 
above, so that in fact we just need to remove the t argument everywhere above because 
all the oscillating phase factors accrued actually cancel. Therefore the final result is: 

dp{t) 



dt 



if) ' 



r 



-ffTLF 



sq-TLF 



Pit)] 



j 



si^'^pmr - pit) 



U) 



s'l>p{t)sf - (s 



pit) + Pit)s 



+ ^- 



)/2 



+ r 



This master equation is valid when flj 3> h'j and f2j ^ maxjA^}, Vj. The first 



inequality is needed in the interaction picture (A.l) and the second inequality is the 



standard requirement for the Born-Markov approximation. 

Within the same framework we can analyze other situations, such as interacting 
TLFs, and additional qubits coupled to the TLFs. The master equation for interacting 
TLFs is the same as above, but with the additional Hamiltonian Vtlf-tlf 

and the requirement that minjrij, fi/j} 3> fJ-j.k- Additional qubits 



can be included under similar conditions. 
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